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ABSTRACT 

We present near-UV transmission spectroscopy of the highly irradiated transiting 
exoplanet WASP- 12b, obtained with the Cosmic Origins Spectrograph (COS) on the 
Hubble Space Telescope (HST). The spectra cover three distinct wavelength ranges: 
NUVA (2539-2580 A); NUVB (2655-2696 A); and NUVC (2770-2811 A). Three in- 
dependent methods all reveal enhanced transit depths attributable to absorption by 
resonance lines of metals in the exosphere of WASP- 12b. Light curves of total counts 
in the NUVA and NUVC wavelength ranges show a detection at a 2.5cr level. We 
detect extra absorption in the Mg II AA2S00 resonance line cores at the 2.8cr level. 
The NUVA, NUVB and NUVC light curves imply effective radii of 2.69+0.24 R^, 
2.18+0.18 R/, and 2.66+0.22 R/ respectively, suggesting the planet is surrounded by 
an absorbing cloud which overfills the Roche lobe. We detect enhanced transit depths 
at the wavelengths of resonance lines of neutral sodium, tin and manganese, and at 
singly ionised ytterbium, scandium, manganese, aluminum, vanadium and magne- 
sium. We also find the statistically expected number of anomalous transit depths at 
wavelengths not associated with any known resonance line. Our data are limited by 
photon noise, but taken as a whole the results are strong evidence for an extended ab- 
sorbing exosphere surrounding the planet. The NUVA data exhibits an early ingress, 
contrary to model expectations; we speculate this could be due to the presence of a 
disk of previously stripped material. 

Subject headings: stars: individual (WASP- 12) 
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1. Introduction 



Observations of the transiting extrasolar planets HD209458b and HD 189733b revealed an 



enhanced transit depth at the wavelengths of several UV resonance lines (Vidal-Madjar et al. 



2003 [ 2004, Lecavelier des Etangs et al.|2010 1. These UV lines from the ground state are sensitive 
probes of the presence of atomic and ionic species. Their presence enhanced the effective radius 
of the planet during transit, implying the planet is surrounded by an extended cloud of size 
comparable to or larger than its Roche lobe ( ] Vidal-Madjar et al.|[2003l |2004| |Ben-Jaffel|[2007 , 



Vidal-Madjar et al.|2008 ). This was attributed to a hydrodynamic 'blow-off' of the planet's outer 
atmosphere caused by the intense irradiation suffered by this hot Jupiter exoplanet. An alternative 
explanation in which the planet is surrounded by a cloud of energetic neutral atoms caused by 
interactions with the host star's stellar wind has, however, been suggested ( Holmstrom et al.|2009 ; 
Ekenback et al.||2010 ). WASP- 12b is one of the hottest and most irradiated transiting exoplanets 
and orbits extremely close to a late F-type host star ( |Hebb et al.|2009 ). WASP- 12b is, therefore, an 
attractive target to explore the properties of the phenomenon observed in HD209458b, and might 
yield evidence distinguishing between the suggested underlying causes. 

The initial UV observations of HD209458b were in the far UV around the Ly a emission line. 
The abundance of hydrogen makes this an attractive line to observe, but the temporal and spatial 
variability of stellar Ly a emission is a highly undesirable complicating factor. For this reason, 
and to obtain better signal to noise, we observed WASP- 12 in the near-UV where there are many 
other resonance lines ( Morton|1991 2000), including the very strong Mg II UV resonance lines. 
This work became possible with the installation of the Cosmic Origins Spectrograph (COS) on the 
Hubble Space Telescope (HST) reinstating and enhancing our capabilities for UV spectroscopy. 



2. Observations and data reduction 

The planet-hosting star WASP- 12 was observed for five consecutive HST orbits on 2009 



September 24th and 25th with COS; see Green et al. (2003), Green et al. (2010, in preparation) 



and Osterman et al. (2010, in preparation) for details of COS. We used the NUV G285M grating 
at the 2676 A setting, which provides non-contiguous spectra over three wavelength ranges 
(NUVA: 2539-2580 A, NUVB: 2655-2696 A, and NUVC: 2770-2811 A) at a spectral resolution 
of i? ~20000, in TIME-TAG mode. The exposure time was 2334 sec in the first HST orbit and 
about 3000 sec per subsequent HST orbit. The optical ephemeris gives ingress during the second 
HST orbit and egress in the fourth HST orbit. 

We downloaded data from M AS iQ adopting CALCOS V.2.11b0for calibration. Despite the 
early date of our observations, the CALCOS reference files used were at a fairly mature stage 
for the NUV data. In particular, the flat field had been updated to its flight version. In our time 
series analysis we used the count rates obtained after background subtraction, rather than the 
flux calibrated spectra. The high quality flat-field and the relatively low background of the NUV 
channel, mean the uncertainties are dominated by poisson statistics. The count rates summed over 
wavelength are roughly 10 counts"'; 28 count s'; and 13 counts"' respectively for the NUVA, 
NUVB, and NUVC ranges. The resulting signal to noise ratio (SNR) per pixel in the NUVB 
spectrum is ~10 for each 3000 sec exposure. 

Figure [T] shows the total summed spectrum in comparison with synthetic fluxes from the 



LLmodels stellar model atmosphere code (Shulyak et al. 2004), assuming the fundamental 



'http : //archive . stsci . edu/ 



^See the COS Data Handbook for more information on CALCOS: 



http : //www . stsci . edu/hst/cos/documents/handbooks/ 



ldatahandbook/COS_longdhbcover . html 



parameters and metallicity given by Hebb et al. (2009). We used the VALD database (Piskunov 
eraL]p95| [Kupka et al.|p99t [Ryabchikova et al.|p99l ) for atomic line parameters and SYNTH3 
( |Kochukhov|2007| ) for spectral synthesis. All three regions are strongly affected by many blended 
photospheric absorption lines; we observe no unabsorbed stellar continuum. The NUVB region 
is closest to the continuum, while the NUVC region is strongly absorbed by the Mg II doublet at 
2795.5 A and 2802.7 A. 



-10- 

3. Detection of a wavelength dependent planet transit 

We expect the planet's atmosphere to absorb particularly in the resonance lines of abundant 
elements. We used three methods to examine the data for wavelength-dependence of the transit 
light curve. 

3.1. The Mg II lines 

The most prominent observed lines in the stellar photosphere are the Mg II lines, and we 
might expect these strong lines to be detectable in the planet's atmosphere too. We adopted the 



method pioneered by Charbonneau et al. (2002) in their detection of the sodium D lines in the 



atmosphere of HD209458b. We divided the NUVC data, which is centred on the Mg II resonance 
hnes, into "blue" (b), "red" (r), and "center" (c) spectral regions. We tried three different widths 
of the center band, "narrow" (n), "medium" (m), and "wide" (w); see Fig. [T]and Table[T} For each 
of these bands we produced a photometric time series, and the associated uncertainty based on 
Poisson statistics. Each photometric index was obtained by summing the observed counts over the 
given wavelength range. In this way "nfo(0" indicates the count rate in the blue side "narrow" set 
at the time t. 
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Fig. 1. — Comparison between the observed mean spectrum of WASP- 12 (thick black line) and 
LLmodels synthetic fluxes (thin red line). The blue dashed line shows the modeled level of the 
stellar continuum flux. The three observed spectral ranges are defined as NUVA, NUVB and 
NUVC from top to bottom. In the bottom panel the vertical lines show the limits applied for 
the wavelength regions, in laboratory wavelengths, accounted to produce the photometric indexes 



described in Sect. 3.1 showing as example the wavelength regions of n^, n^ , and n^. 
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The stellar limb darkening could potentially cause a color-dependent transit shape (e.g, 
Brown et al.|2001 1. To assess this we calculated the difference of the blue and red spectral regions 



for the "n", "m", and "w" bands as a function of time (see Eq. 1 of Charbonneau et al^|2002 ). 



We looked for variations in the transit depth due to the stellar limb darkening calculating the 
difference between the mean photometric indexes obtained in- and out-of-transit (see Eq. 2 of 
Charbonneau et al.|2002 ). All values we obtained were clearly consistent with no variation. 



To examine time dependence using Charbonneau et al. (2002 1 's method we calculated in 



each band ("n", "m", and "w") the difference between the light curve of the central band and the 
mean light curve of the blue and red bands: 

riMgit) = ndt) - [Uhit) + Urit)] I 2 

T^Mgit) = tndt) - [mh(t) + nirit)] / 2 (1) 

WMg{t) = Wc(t) - {Whit) + Writ)] I 2. 

In this way, we removed any limb darkening dependence. Again, the time series have 
RMS scatter consistent with photon noise: {(T\nMg{tout)] ~ o-\mMg{tout)\ ~ o-\w Mgitout)] ~ 
3.4 X 10"^ counts"^). We then calculated the difference between the mean in-transit and 

Table 1 : Limits adopted to define the analysed wavelength regions around the Mg II resonance 

lines. 

Band Wavelength Band Wavelength Band Wavelength 

range [A] range [A] range [A] 

Hb 2782.75-2795 m^ 2782.75-2787.5 Wh 2782.75-2783.75 

rir 2805-2817.25 m, 2812.5-2817.25 w, 2816.25-2817.25 

ric 2795-2805 m, 2787.5-2812.5 w, 2783.75-2816.25 
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out-of-transit flux: 



^fiMg = riMgiUn) - riMgitout) = (3.5 + 4.1) X 10 ''counts 1 

AmMg = mMgitin) - niMgitout) = (-4.7 + 4.1) x 10"^ counts"' (2) 



AWMg = WMgiUn) - WMgitout) = (-1 1 .4 + 4.1) X 10 ^ COUUt S ^ 

These results show the detection of a deeper transit in the "m" and "w" bands at l.lcr and 2.8cr, 
respectively. Since the value obtained in the "n" band is comparable to the resulting photon noise 
error bar we believe that the non-detection is due to the very low signal level in n^ along with 
absorption occurring in the wide n,- and nj, bands. The size and the significance of the detection 
increases as the signal included in the center band increases, just as we would expect if the 
enhanced transit depth in the Mg II doublet is genuine. 

3.2. The transit light curve 

We compared the light curves obtained for each observed wavelength range and the one 
calculated from visible photometry, as shown in Fig. [2} 
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The NUVB wavelength range is the closest to the continuum and shows a transit depth that 
matches, at ~ Icr, the transit light curve derived by Hebb et al. ( |2009 1 from optical photometry. In 



the NUVA and NUVC wavelength ranges we obtained a deeper transit at about 2.5cr level. These 
three light curves were normalised to the line passing through the out-of-transit photometric points 
(first and fifth exposures). The slope of the three normalisation lines are 3.8 x 10"^ for the NUVA 
region, 3.3 x 10~^ for the NUVB region, and 1.0 x 10'^ for the NUVC region. These values are 
small enough that the applied normalisation did not change the transit shape. 

The NUVC spectral region is clearly dominated by the Mg II resonance lines that are likely 
to be responsible for the detected extra depth in the transit light curve. The NUVA spectral region 
includes resonance lines of Na I, Al I, Sc II, Mn II, Fe I, and Co I ( |Morton|[T99T| [2000). The 



stellar spectrum is dominated by Mg I and Fe I lines coming from low energy levels. Probably, 
these spectral features, likely to be present also in the spectrum of the planet atmosphere, produce 
the observed deeper transit (see |Vidal-Madjar et al.|2004 , for a similar case). 



The end of the second exposure is at the phase of the planet ingress, as shown in Figure [2} It 
is notable that the NUVA flux during the second exposure lies below the out-of-transit level by 
~ 2cr. We divided this particular exposure into three equal sub-exposures plotted as black crosses. 
These suggest an early ingress in the NUVA spectral region. 

3.3. Detection of other elements 

In each of the three observed wavelength ranges we calculated a ratio spectrum (cIa) between 
the in-transit spectrum (m^) measured in the third exposure and the out-of-transit spectrum (out;i), 
the mean of the first and fifth exposures. To these ratio spectra we associated two different 
uncertainties: (i) the standard deviation from the mean, d, which we denote cr^/j|^„. (ii) The 
uncertainty for each individual wavelength point in the ratio spectrum from the propagated 
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Fig. 2. — Light curve obtained for each observed wavelength range (NUVA: open black circles 
- NUVB: open red squares - NUVC: open blue triangles). The horizontal error bar defines the 
orbital phase range covered by each observation. The vertical uncertainty comes from a Poissonian 
treatment of the error bars. The full green line shows the MCMC fit to the optical transit light 
curve ( |Hebb et al.|2009| ). The black crosses show the NUVA spectral range split into three equally 
exposed sub-exposures. Lines indicate the normalisation gradient applied. 
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uncertainties. We denote this cr^,|j„^. Expressed symbolically: 



dx = —— (3) 

out^ 



and 

where // is the number of points, (Tin^ = V"^: and cr„„,, = ^^ouf^. In NUVA, NUVB, and NUVC 
(TdA\e:,p is 0.34, 0.12, and 0.76 respectively. crdA\,,rop varies with wavelength, as shown in Fig. [3} 

Table |2] lists the wavelength points of J^ (in laboratory wavelengths) with deviations of more 
than 3cr from d, assuming both a = crj,|„ (left column) and a = cr^^i ,„, (right column). Assuming 
a Gaussian distribution and having N= 1024x3 wavelength points, we expect 9 points in the d^ 
array to fall outside 3cr from the mean. Since the number of detected deviating wavelength points 
is much larger than nine we looked for correspondences with resonance lines ( Morton|[1991 



2000). Table |2] lists the deviating wavelength points and the corresponding resonance lines. We 



include occurrences of resonance lines within a few km s ' of a deviating wavelength point, for 
example the Sc II line at 2563.190 A. 
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In the NUVA wavelength region and adopting cr = cr^^i^^^ we obtained 3cr deviations 
corresponding to the position of three resonance lines: Yb II at 2538.662 A, Sc II at 2540.822 A, 
and Mn II at 2576. 106 A. Assuming instead cr = CTd.i^,,,, Sc II at 2540.822 A and the Na II doublet 
at 2543.8 A are picked out. In the NUVB region we only find the V II line at 2683.090 A and 
only assuming cr = crj^i^,- However three other V II resonance lines and an Al II line lie close 
to other detected deviating points. In the NUVC region we recognize immediately that most of 
the deviating points are in the core of the Mg II resonance lines, both assuming cr = crd^,\^„ and 
cr = crrf,|,„„p. We also pick out the Mn I line at 2801.082 A, while the Mn I line at 2798.269 A lies 
close to the wavelength of another deviating point. 

Figure [3] shows the cores of the Mg II resonance lines. We show the observed spectrum, 
d^, cr^ii,,,,,,,, and the deviating wavelength points both assuming cr = cr^/^,^ and cr = cr^/,|^^^_^^. With 
cr = crdj^i^ the deviating points correspond to the core of the Mg II line where the signal level is 
low. This is to be expected: the low count rates at these wavelengths lead dA to be very noisy here. 
In contrast, with cr = crii;i\„.„, each element of the dx spectrum is assessed against its own Poisson 
error. In this case the deviating points are all below the mean rate spectrum, and the deviating 
points appear at the margins of the line core. These points indicate excess Mg II absorption during 
transit. This is attributable to absorption by the planet's atmosphere. This pattern is seen not only 
for the two Mg II resonance lines, but also for the Sc II line at 2563.190 A. This line, together with 
the Mn I line at 2798.269 A, has the intriguing property that the difference between the position of 
the resonance line and of the detected deviating wavelength point(s) corresponds to a velocity of 



-30 km s ' (about 3 resolution elements), close to the planet escape velocity of ~37 km s ^ (Hebb 



et al.||2009 ), although it would not then be clear why this pattern does not appear also for other 



detected lines of the same ion. 
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Table 2: Wavelength of the spectral points deviating more than 3cr, adopting two different crs: 
cTrf^l^^p and cTrf^ip^p. For each detected deviating point we show the resonance line found lying at 
the same position or close to it (*). In the NUVB region we did not detect any deviating point 
assuming cr = cr^^|,„,. The deviating points marked with a # deviate by >3.5cr from the mean. 
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There are 95 known resonance lines lying within the observed wavelength ranges, including 
those of heavy elements. This is a small minority of the more than 4000 lines present in the stellar 
spectrum. The fact that we find deviating points predominantly at wavelengths corresponding to 
resonance lines strongly suggests we are detecting features produced by the planet atmosphere. 
Reassuringly, with either definition of cr, we obtained nine points that do not match any known 
resonance line, in perfect accordance with statistical expectations. We repeated the exercise 
picking out deviations in excess of 3.5cr, obtaining almost the same deviating points at the position 
of known resonance lines and fewer points where no resonance lines were found. 

4. Discussion 

We have performed three independent analyses, each of which suggests absorption in 
the resonance lines of metals from an extended atmosphere surrounding the transiting planet 



WASP- 12b. In Section 3.1 we found a deeper transit in the core region of the Mg II doublet at the 
2.8cr level. 



In Sect. |X2} the transit depths in the NUVA, NUVB, and NUVC wavelength ranges 
respectively imply eff'ective planet radii of 2.69+0.24 Ry, 2.18+0.18 R/, and 2.66+0.22R^. 
WASP-12b's optical radius is Rp = 1.79 + 0.09 Rj while the mean Roche lobe radius is 2.36 Ry 
using Paczyhski's (1971) prescription. 

Table |2] shows that we detect enhanced transit depths at the wavelengths of resonance lines 
of neutral sodium, tin and manganese, and at singly ionised ytterbium, scandium, manganese, 
aluminum, vanadium and magnesium. Finally we detect an enhanced transit depth within 0.1 2A 
of a resonance line of doubly ionised europium. We also find the statistically expected number of 
anomalous transit depths at wavelengths not associated with any known resonance line. 

Taken as a whole, these results constitute compelling evidence that WASP- 12b is surrounded 
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Fig. 3. — The black line shows the observed spectrum obtained averaging the five available COS 
spectra. The red line shows the d^ spectrum, magnified five times and shifted upwards for display 
reasons. The blue lines show the cr^^i ^,.^ , spectrum (full line) and the values of cr^^i^, (dashed lines). 
The full black circles show the position of the deviating points assuming cr = cr^^|„,, while the full 
green triangles show the position of the deviating points assuming cr = (Td^\ „ . 
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by an exosphere which over-fills the planet's Roche lobe, confirming predictions by |Li et al. 



( |2010 ). This exosphere is likely composed of a number of elements/ions, including probably 



Na I, Mg I, Mg II, Al I, Sc II, Mn II, Fe I, and Co I. The phenomenon found in HD209458b 
( |Vidal-Madjar et aL]|2003[ |2008[ ) probably occurs generally for hot Jupiter exoplanets. By 



analogy with HD209458b, and as WASP- 12b and its host star are almost certainly predominantly 
composed of hydrogen, we expect that this exosphere is hydrogen rich. 



Models by Yelle (2004) suggest that elements other then H and He should not be present 
in the upper atmosphere due to the low vertical mixing rate, but this takes Jupiter as the starting 
point. WASP- 12b is extremely close to the host star and consequently the stellar irradiation and 
tidal effects could induce prodigious mixing, affecting the chemistry of the planet atmosphere. 
Our detections of several metallic elements and/or ions is certainly consistent with a metal-rich 
atmosphere for WASP- 12b. 

The most surprising result is provided by the juxtaposition of our data with the optical 



ephemeris. We took contemporaneous optical photometry with OU-OAM PIRATE (Kolb et al 



2009| ) which showed the ephemeris of |Hebb et al. (2009) remains accurate. Figure [2] shows the 



NUVA transit has an early ingress and an egress consistent with the optical ephemeris. In contrast, 
naive momentum considerations and hydrodynamic simulations would instead suggest that the 
effect of a diffuse cloud surrounding the planet would be to smear and delay egress while ingress 



is relatively unaffected, see e.g. Fig. 1 and 2 of Schneiter et al. (2007 ). 



In detail, the shape of the diffuse cloud may well be element/ion dependent since different 
elements/ions behave differently in the presence of strong radiation pressure. This can explain why 
we observe different transit shapes in the NUVA region and the other regions. As Fig. [T] shows, 
the stellar spectrum in the NUVA region is strongly absorbed by a plethora of lines, dominated 
those of neutral elements. The NUVC region is also strongly absorbed in the stellar photosphere 
but predominantly from the Mg II doublet. It is presumably the cumulative absorption from many 
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relatively weak spectral lines in the planet's exosphere which creates the excess transit depth in 
the NUVA region, while Table [2] and Eq. |2] demonstrate that planet's absorption in the NUVC 
region is associated with the Mg II doublet. The Mg II ion will experience different forces to 
neutral atoms in an environment where there is certainly a strong radiation field, and strong and 
varying large-scale magnetic fields are also likely. The NUVB light curve is least deviant from 
the optical transit, and this is consistent with the relative dearth of strongly absorbing lines in this 
spectral window, c.f. Fig.[T] 

We do not have any detailed explanation for the observed early ingress in NUVA, but we 
speculate the effect could be produced if material is lost from the planet exosphere and forms a 
diffuse ring or torus around the star enveloping the planet's orbital path, as models suggest ( |Li | 



et al.||2010 ). The orbital motion of the planet through this medium might compress the material 
in front of it. This could increase the opacity of the medium through which the star is viewed 
immediately before first contact. A void in the medium might be expected to form behind the 
planet, and consequently the egress is relatively unaffected by the diffuse ring. 

Our observations demonstrate that COS spectroscopy of transiting exoplanets has the 
potential to detect many species via transmission spectroscopy, and to measure velocities and 
deduce spatial distributions. There are now about 40 known transiting exoplanets with orbital 
periods shorter than that of HD209458b. Many of these transit stars significantly brighter 
than WASP- 12b. COS spectroscopy of brighter examples will allow us probe the exosphere 
species-by-species examining their density, velocity and spatial distributions. This detailed 
information should allow us to determine whether these planets really are being photo-evaporated 
by their host stars, and, if so, to empirically deduce the mass loss rate. We encourage detailed 
element/ion dependent modeling of the exosphere in the highly irradiated environment of 
WASP-12b and similar systems, and observations of other similar extrasolar planets. There is a 
rich new parameter space to explore! 
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